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1 Problem with a large initial gradient 


A simplest model of the motion of continuum, which takes into account nonlinear effects and 
dissipation, is the equation of nonlinear diffusion 


du du d‘^u 
dt ^ dx ^ dx"^ ’ 


£ > 0 , 


( 1 . 1 ) 


for the hrst time presented by J. Burgers [T]. This equation is used in studying the evolution 
of a wide class of physical systems and probabilistic processes, for example, acoustic waves 
in fluid and gas [2]. 

Let us consider the Cauchy problem 


■u(x, 0, £, p) = A(a;p ^), f = 0, a; G M, p > 0, 


( 1 . 2 ) 


where p is another small parameter. The initial function A is smooth and it has hnite limits 

Aq = lim A(s). 

We see that the origin (a: = 0 and t = 0) must be a singular point in the problem 
under consideration. It is known that the behavior of solutions of differential equations with 
a small parameter at a higher derivative in a neighborhood of a singular point sometimes 
becomes self-similar. Then it is effective to analyze the solution using the renormalization 
group method [3]. The success of applying renormalization is conhrmed by examples from 
hydrodynamics |1], mechanics [5], and asymptotic analysis of solutions of nonlinear parabolic 
equations [6]. 

In paper |^, for the solution of a more general quasi-linear parabolic equation 


du d(p{u) 
dt dx 


= e 


di^u 

dx"^' 


t ^ 0, 


1 





with condition fll.2jl the following asymptotic formula is obtained: 


u(x,t,E,p) 



(1.3) 


as £ —)• 0 and p/e —?■ 0 uniformly in the strip {{x,t) : a: G M, 0 ^ t ^ T}, where T is the 
solution of the limit problem in the inner variables rj = x/e, 9 = t/e 


ar ay(r) 

do dr] 


d^T 

drj'^ 


0, r(p,o) 


Aq , V <0, 
A+, p>0. 


The procedure of obtaining expression fll.3|) can be also seen below in details on the example 
of the KdV equation. The difference consists in the scales of the inner variables. For the 
solution of the Burgers equation fll.ip as £ —)■ 0 and p/e —)■ 0 in the strip {(x, f) : x G M, 0 ^ 
t ^T} there holds the asymptotic formula 


u{x,t,e,p) = 


Ao"' - Ao 


Aq exp 


t(A+)^ - 2A+X A+ps 
4e 2e 


erfc 


A+f 


X + ps 


2y/ei 


+ 


Aq exp 


^(Aq) -2AoX Agps 


4e 


+ 


2e 


erfc 


X 


exp 


(( A +)2 _ 2 A +_£ ^ ^ I 


4e 


2e 


X — Aq f — ps 
2\/ef 

Aq t — X + ps 


X 


2y/ei 


+ 


exp 


^(Aq) - 2 AoX AqPs 


4e 


+ 


2e 


erfc 


X — Aq f — ps 
2^/ei 


-\ -1 


ds + O ■ 


2 Renormalized weak discontinuity 


In the present paper, we give an asymptotic solution of the Burgers equation 


du du di^u 

dt dx dx"^ 

with the initial condition 

m(x, 0) = eA j f = 0, X G M, 

where A(s) —)■ — s as s —)■ —oo and A(s) —)■ 0 as s —)■ +cxd. 

By analogy with fll.3p we postulate the expression 


R{x, t, e) 


OO 

j A"[SWX-EE, t) is, 


— OO 


( 2 . 1 ) 


( 2 . 2 ) 
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where uq = —2^x/'^ is the solution of the limit problem with a weak discontinuity |8], 


duo duo d'^uo \-x, X < 0, 

+ ““^=9^- =i 0. .SO, 


dt 


d' is the solution of the heat equation = ^xx- Thus, we obtain a renormalized asymptotic 
solution in the form 


R{x,t,e) = — / A"'(s) ln\['(x — £s, t) ds. 


Theorem 1. Let A{s) be a thrice differentiable function and let 


Then the function 


- l + 0(|s|-3), CX), 

(2.3) 

‘‘ff = 0(s-=), s ^ +00, 

ds 

(2.4) 

d‘^A{s) _4 

ds^ =0(. ), .^oo. 

(2.6) 

OO 

^ 2 /• d^A(s) , , , 

R(x, t,e) = — / —-^-2- In T(a; — es, t) ds 
s J ds'^ 



is an asymptotic solution of the Burgers equation fl2.ip . where 


'^{x,t) = / exp 


(x - a)^ 




Proof. Substituting u = R{x,t,e) into fl2.ip . we hnd 

CXD CXD 

Ri + RR^ - =-f j h"'{s) ds J A"'{s)\n'^{x - Es,t)ds- 


_2 f 

£ J '^'^{x — es,t) 

— OO 

Let us esimate the integrals on the right-hand side. From fl2.3p - fl2.5l) it follows that 


l/Ve 1/Ve 

J A"'{s)ds = 0{e^), J sA'"{s)ds =—1 + 
-l/^/£ -l/Ve 
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Then we obtain 


l/v/e 

A'"(s) ln\I/(a; — £s,t) (is = f A'"(s) 

-l/v/i 


lnT(a;,t) — £s 


t) 


ds + 0{e'^) = 




T(a;, t) 


, ,\ 1 /3.(37 — ES.t) I ^ XX 

A"(s)^ - -A ds = ' 


^(x — es, t) 


l/Vs 

111 [ sA"'{s)ds + 0{e^) = 




-i/Vi 


= e 








— es, t) 

— 00 

= 2e 

Consequently, we conclude that 


Vl/2 ^3 


l/\/e 


sA'''{s)ds + 0(£^) = 


-l/^/e 


\]>3 

^ XX ^ X 


Vl>2 ^3 


+ 0(£^). 


Rt + RRx — Rxx — 0(e). 

Since R{x,0,e) = eA(x/e), by the continuity of the inverse operator, we can also state 
that R{x, t, e) is an asymptotic approximation of the exact solution of the Cauchy problem 

(EiD-dii. 


3 Korteweg—de Vries equation 


Studying properties of solutions of the Korteweg-de Vries equation 


du du d^u 


t ^ 0, £ > 0, 


(3.1) 


is of indisputable interest for describing nonlinear wave phenomena |2]. Open questions con¬ 
cerning the behavior of solutions are still the subject of attention for modern researches [9]. 
Some mathematical results about the solution of the problem in various cases can be found 
in 110] and [TT] . 

Consider the Cauchy problem for fld.ip with the initial condition 


u{x,Q,e, p) = A[xp ^), t = 0, a; G M, p > 0. (3.2) 

Here, we assume that the initial function A is bounded and possesses the hrst derivative, 
which sufficiently fast tends to zero at inhnity. It is required to hnd an asymptotic approx¬ 
imation of the solution u{x, t, e, p) to problem fl3.ip - fl3.2p in small parameters e and p for 
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all X and finite valnes of t. It is clear that the structnre of the asymptotics must depend on 
the relation between parameters £ and p. We assume the fulhllment of the condition 

p = ^ 0. 

In paper na, a uniformly suitable asymptotic approximation of the solution to problem 
fl3.ip - fl3.2p is constructed using the renormalization method in the following most simple 
form. Let us pass to the inner variables 

X = \/i?7, t = y/e9, 


since this allows one to take into account all terms in equation fl3.ip . As a hrst approximation, 
we take the solution of the equation 


dZ dZ d^Z 

d9 dr] drj^ 


(3.3) 


with the initial condition 


Z{r],0) 


Aq, V<0, 
Aq , h > 0, 


(3.4) 


where A^ = lim A(s). As a model of collisionless shock waves, problem fl3.3^ - 113.4p was 

s^zLoo 

studied by A.V. Gurevich and L.P. Pitaevskii in [T3] . 

Let us construct the expansion of the solution in the following form: 


u{x,t,e, p) = Z{r],9) + pW{r],9, p) + a > 0, (3.5) 

where the addend pW{p, 9, p) must eliminate the singularity of Z at the initial moment of 
time. Then the function W satishes the linear equation 


dW diZW) d^W ^ 

d9 dp dp^ 

Differentiating equation fl3.3p with respect to p, we hnd that the expression 


(3.6) 


G{p,9) = 


dZ{p,9) 




An 


dp 


satishes equation fl3.6p . Moreover, G is the Green function, because 


OO OO 

^ - ^+ J Z{p, 0)fip) dp = /(O) 

— OO —OO 


for any smooth function / with compact support, thus G{p, 0) = S{p). 

Let us choose the solution W in the form the convolution with the Green function G 
so that the asymptotic approximation would satisfy the initial condition fl3.2p . As a result, 
expansion fl3.5p becomes 


u{x, t, e, p) = Uo{x, t, e, p) + 0(/i“), 
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where 


Uo{x,t,e,p) = Z{r],9) + 


/i 


- Ao" 


dZ{r] — ps, 9) 
dr] 


[A(s) — Z{s, 0)] ds. 


Integrating by parts, we obtain the asymptotic approximation 


u{x,t,e,p) Uo{x,t,e,p) = 


A^-Ao 


X — ps t 


A'(s) ds. 


A rigorous mathematical justihcation of this formal representation is out of the framework 
of the present paper. It is clear that the behavior of the convolution integral is essentially 
determined by the solution of fl3.3p - fl3.4p with discontinuous initial data that conhrms the 
importance of the Gurevich-Pitaevskii problem. 

Constructing complete asymptotic expansions of the solution near the singular point by 
the standard matching method may be connected with serious difficulties. In fact, it is neces¬ 
sary to solve the scattering problem for a recurrence system of partial differential equations 
with variable coefficients [2]. In addition, the investigation of the shock wave generated 
by gradient catastrophe shows that the asymptotics of the solution in a neighborhood of a 
singular point may have a multiscale structure |8]. 

The renormalization approach allows one to construct a uniformly suitable asymptotics 
in the whole domain of independent variables avoiding difficulties arising from the matching 
procedure. 

In particular, for A(|" = 0 and Aq = a > 0 the following formula was obtained in 
paper na: 

/ N /X + at\ . /'x — 2at/3\ 

u{x,t,e,p) ^2Ai - 1-A' ' ' 


P 


P 


■) 


at 

P 


2/3 


A' 


-1 


X — aty 
P 


2dffi 






dy, 


where dn(r;, a) is the elliptic Jacobi function 
dn(M, m) = sjl 


ip(u) 


dv 


msm (p. 


u = 


<^{y) = 




//--[! + a^{y)] \ , 1 + a 


J \/l — m sin^ V 
0 

2(T^(1 — CT^)iP(cr) 
E{a) — (1 — a‘^)K{a) 


= ^y, 


K{a) and E{a) are complete elliptic integrals of hrst and second kind: 

7r/2 7r/2 

E{a) = / \/l — sin^ v dv. 


= j 


dv 


\/l — cr^ sin^ V 


0 ’ 0 

Note that using the change x = p^, t = pr for fl3.ip - fl3.2l) . we obtain the problem 


du du 


-2 


d^u 


dr^^'di^^ de 


= 0 , u\r=o = HO, 


which formally contains only one parameter. However, in this case one should study asymp¬ 
totics as /i“^ —)■ oo and r —)■ cxd. 
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4 Perturbation of a weak discontinuity 


Now consider the Cauchy problem for the KdV equation 


du du 


with the initial condition 


d^u 

dx^ 


= 0 


u{x^ 0, e) = eA j 


t = 0, X G 


(4.1) 


(4.2) 


where A(s) = —s0(—s) + 0(s ^) as s ^ oo, 0 is the Heaviside function. By analogy with 
the case of the Burgers equation we postulate the expression 


R{x, t, e) 


OO 



— CO 


es, t) ds 


(4.3) 


as an asymptotic solution for 0 < f ^ 5 < 1, where $(x, f) is a smooth solution of the 
Faminskii problem [15] 


= 0, <h(x, 0) = -x0(-x) (4.4) 

with an initial weak discontinuity also called a contact discontinuity. It is easy to show that 
the function fl4.3l) exactly satishes the initial condition fl4.2p . 


Theorem 2. Let A(s) be a twice differentiable function and let 

A(s) = —s 0(—s) + 0(s“^), = —0(—g) + 0(|s|~^), s ^ OO. (4.5) 

ds 

Then the function fl4.3p is an asymptotic solution of the Korteweg-de Vries equation fl4.ip . 
Proof. From fl4.5p it follows that 


l/Vi 

1/v/e 

/ A" {s)ds = 1 + 

/ sA"{s)ds = 0{e) 

l/Vi 

-i/V^ 


Taking into account these relations, using fl4.4p . and substituting u = R{x,t,e) into fl4.ip . 
we hnd 


Rt + RRx + Rxxx = / A"{s)^{x — es,t) ds / A"{s)^x{x — es,t) ds— 


OO 

j A"(s)<h(x — £s, f)^x{x — £s, t) ds = 

— OO 
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